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Abstract 

The preliminary data on hadron multiplicities measured in central Pb+Pb 
collisions at 158Y GeV/c are analyzed. The ideal hadron gas model fails 
to give a reasonable explanation to the Pb+Pb data sets. We study the 
possible effects of pion enhancement due to different hard-core repulsion for 
pions and other hadrons and strangeness suppression because of incomplete 
chemical equilibrium. Each of these two modifications improves the results. 
The combined effect of these two mechanisms leads to an extremely good 
agreement with the data. An interpretation of the obtained results in terms 
of the possible quark-gluon plasma formation at the early stage of the collision 
is also discussed. 
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I. INTRODUCTION 


Different versions of thermal hadron gas (HG) models have recently been used to £t 
the data of particle number ratios in high energy nucleus-nucleus {A+A) collisions at BNL 
AGS and GERN SPS (see e.g. [1-8]). These particle number ratios were measured in fixed 
rapidity intervals. Particle abundances in HG models are described in terms of two chemical 
freeze-out parameters, temperature T and baryonic chemical potential hb- However, the 
restriction on hadron yields in a fixed rapidity interval causes some problems. The universal 
energy distributions of all hadrons in HG models are transformed into different rapidity 
distributions for different hadron species — the width of rapidity distribution is decreased 
with increasing particle mass. Therefore, to calculate these particle number ratios measured 
in the fixed rapidity interval, additional model assumptions on collective flow effects and 
specific values of thermal freeze-out parameters are required. 

In the present work, we analyze the data on hadron multiplicities in central Pb-|-Pb 
collisions at 158A GeV/c measured at GERN SPS. The aim of our thermal HG analysis of 
hadron multiplicities in a full 47r-geometry is to avoid the above-mentioned complications. 
The basic assumption here is rather strong that the values of the chemical freeze-out pa¬ 
rameters T and Pb in different elements of the system are equal. However, the collective 
flow effects are still permitted and the validity of a pure thermal (“fireball”) model is not 
necessarily required. 


II. THE MODELS 


The ideal gas equation of state for a mixture of different particle species i=l,...,h with 
chemical potentials p.i,...,/i/i is given by the pressure function in the grand canonical ensemble 
as 
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where the integration variable k is the momentum, and di and rrii are the degeneracy factor 
and the mass of the zth hadron, respectively. The distribution function is written as 
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The statistic factor rji is —1 for bosons and +1 for fermions. The thermal number of ith 
particle species is expressed as 

= y^ ^ ^ j'dkeMk), (3) 

in which V is the total volume of the hadron system assumed to be equal to the sum of 
the proper volume elements at the chemical freezeout. Hadron chemical potentials /ij are 
expressed as 


/ii = bifiB + Sifis (4) 

in terms of baryonic {fis) and strange {fis) chemical potentials with hadron charges bi = 
0, ±1 and Si = 0, =f1, =f2, t 3. Direct (thermal) particle yields are calculated according to 
Eq. for all known stable hadrons and resonances with mass up to 2 GeV. The total 
multiplicity of ith hadron is the sum of its direct yield and all possible strong decay 
contributions from resonances. This is true for all HG models discussed in the present 
paper. 

There are serious theoretical doubts in the validity of the ideal HG model at the chemical 
freezeout. Particle number densities and total energy density of the HG for chemical freeze¬ 
out parameters found from htting experimental data look artihcially large. This is hardly 
consistent with the picture of a gas of point-like non-interacting hadrons. The introduction 
of the “hard-core” repulsion in the HG model has been widely discussed in the literature 
(e.g., [1-8]). We follow the excluded volume procedure of Ref. 0 with the assumption of 
substituting the total volume V hj V — Yh ViNi. The proper volume of each hadron Vi is 
expressed in terms of its “hard-core” radius ri as Vi = yrf [^]. The HG equation of state 
with Van der Waals (VDW) repulsion can then be written as 
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and particle multiplicities are given as 
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All HG multiplicities (|) are suppressed in comparison with the ideal HG results (j^). The 
same values of Ni in the VDW HG as in the ideal HG correspond to a larger value of volume 
V in the VDW HG than in the ideal HG. An artihcially high energy density in the gas of 
non-interacting hadrons is therefore removed [P,pT[1. 

The ideal HG model can not £t existing data of hadron multiplicities. Numerous attempts 
to £t data on particle number ratios in A+A collisions at AGS and SPS energies have 
revealed two possible ways to modify the ideal HG model: pion enhancement or strangeness 
suppression. In other words, in order to £t data, one usually needs either more pious or 
fewer strange particles in comparison with the ideal HG results. The introduction of the 
phenomenological parameters r* could change particle number ratios in comparison with 
ideal HG results. To enhance pion/hadron ratios we assume a smaller “hard-core” radius 
for pion than those for all other hadrons, < r* = r The alternative dynamical 

explanations of the observed pion excess are presented in Refs. [^,0 in the framework of 
hydrodynamical and transport approaches. 


A chemical non-equilibrium HG with a strangeness suppression was suggested in Ref. [13 


It is usually formulated in terms of a phenomenological parameter < 1, which can be 
connected to additional chemical potential regulating the absolute value of total number of 
strange quarks and antiquarks in the HG system |^. The form of Eqs. (|^,|^) remains the 
same in this approach, but the distribution functions are modihed as 


flik) = 
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where |S'j| is the sum of the number of strange quarks and antiquarks in hadron i. Hence 
Eq. becomes 
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Hadron multiplicity data for CERN SPS nucleus-nucleus collisions were recently fitted with 
Eq. (H) in Ref. 

Hadron multiplicities in HG models depend on 3 thermodynamical parameters, 

V, T, and (strange chemical potential is always dehned by the requirement of zero 
total strangeness). The model with strangeness suppression (H) introduces an additional 
phenomenological parameter jg. Meanwhile, in the VDW HG model (j^), one new parameter 
r is added to the ideal HG formalism — = 0 is assumed for simplicity Note that, 

by putting 7s = 1 and r = 0, the ideal HG model is recovered: Eqs. (H) and are reduced 
to Eq. (I). 

We consider also the HG model which combines both strangeness supperession and hard¬ 
core repulsion effects. The pressure function p in this model is dehned by the equation 



V V y k'ly •••) k'h)y k'i — h'i V y 
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where p® in Eq. (P) is given by Eq. (i, but with distribution functions // (0) instead of fi 
(0). Thermal hadron multiplicities W are then given by 
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III. FITTING PROCEDURE AND RESULTS 

The data for hadron multiplicities in Pb-|-Pb collisions at 158A GeV/c and references to 
the original papers [16-21] are presented in the last two columns of Table I. Because the mea¬ 
surement of hadron multiplicities in Pb-|-Pb collisions at SPS energies are not yet complete, 
we also include several hadron ratios in the central rapidity region, K+/K“, p/p, H/5 and 
U/U (the controversial A/A ratio is ignored) in our htting. Due to equal mass, the kinematic 
complications mentioned in previous paragraphs are not expected to be very important. We 
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assume, therefore, that these ratios in the central rapidity region are approximately equal 
to the corresponding ratios of the total hadron multiplicities. 


The htting procedure is to minimize 

N 

= E 




.=1 i^yrn^ ■ 

To compare the htting of different models to the data, we normalize by 


( 11 ) 


X 

dof 


X 


N — n 


( 12 ) 


where N is the number of experimental data points htted and n is the number of parameters 
in each model. The best ht with ideal HG model is presented in the hrst column (ID) in Table 
I. It has been known that the ideal HG model can not ht all hadron ratios simultaneously. 
This problem has been discovered in the previous attempts to ht data on particle number 
ratios in collisions of S nuclei with diherent targets at GERN SPS p|,^. From our ht of the 
data set in Table I for Pb+Pb it follows that the main problems for the ideal HG model are 
indeed a dehciency of pions, the major contribution to h~, and a surplus of strange particles 
— K° and (j) substantially exceed their experimental estimates. It is possible to improve 
the agreement with h~ and 0 within the ideal HG model (h“=603, 0=6.38) by choosing a 
smaller temperature (T = 125 MeV) and a larger volume {V = 15900 fm^). To keep Nb 
close to data (iV5=370) one needs then a larger value of the baryonic chemical potential 
(ps = 300 MeV). These “small” T and “large” ps parameters in ID HG lead, however, to 
completely wrong antibaryon to baryon ratios: p/p, S/S and D/D all become an order of 
magnitude smaller than their measured values. 

The agreement is improved if one adopts either the strangeness suppression scheme (SS) 
or hadron “hard-core” repulsion scheme (VDW) with = 0 and r > 0 which leads to the 
pion enhancement as mentioned before. The VDW results are shown in the second column of 
Table I. Note that h~ becomes closer to data, but K° and 0 are still too large. The SS results 
are shown in the third column of Table I. The quality of the ht is signihcantly improved. 
However, there is still a dehciency of pions in the SS HG model. Finally, the results of 
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taking both SS and VDW together into consideration are given in the fourth column of 
Table I. The of this £t is remarkably small as shown in the table. In addition, we present 
in Table II all relevant total hadron multiplicities calculated with parameters found from 
SS+VDW. Similar to Table I the hadron multiplicities presented in Table II include both 
direct (thermal) yield (^) and all possible contributions from strong resonanse decays. 


IV. DISCUSSION 


Let us discuss a possible physical interpretation of the obtained results in terms of the 
quark-gluon plasma (QGP) formation at the early stage of A+A collisions. The “enhance¬ 
ment” in strange hadron productions was suggested as a signal for the formation of QGP 
in A+A collisions (see e.g. Ref. One expects that the strangeness equilibration time 

in the QGP is much shorter than in the HG. Therefore, rapid production and equilibration 
of strangeness is expected in the QGP initial state, and the strangeness suppression factor 
( 7 s < 1) is usually associated with an incomplete strangeness equilibrium in the pure hadron 
matter. 

A real picture could be, however, very different. We found that, for the chemical freeze- 
out parameters T and fiB shown in Table I, the strangeness to entropy ratio is larger in the 
equilibrium HG than in the equilibrium QGP. Therefore, 7s < 1 strangeness suppression in 
the HG would become a signal for the formation of QGP at the early stage of A+A collisions 
at GERN SPS energies. The same conclusion was obtained in Ref. . 


Strange quark-antiquark pairs may be primarily produced by hard (with typical momenta 
of several GeV) (anti)quark and gluon collisions at the early non-equilibrium stage of A+A 
at the SPS. If the QGP is formed one can expect the complete strangeness equilibration 
(7s^ = 1 in the QGP). To estimate the strangeness/entropy ratio in the QGP, let us use the 
ideal gas approximation of massless u-, d-(anti)quarks, gluons and strange (anti)quarks with 
rus = 150 MeV. In the wide range of T and fiB {T = 160-300 MeV, /is = 0-400 MeV), one 
buds an almost constant value of the ratio of total number of strange quarks and antiquarks 
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Ns + Ns to the total entropy S: 


Ns + Ns 

Rs = = 0.022-0.025 . (13) 

o 

There is no realistic model for the QGP hadronization. We do not expect, however, 
essential additional contribution to the strange-antistrange hadron production during the 
QGP hadronization: at this stage the typical thermal (anti)quark and gluon momenta are 
only a several hundred of MeV. Therefore, we assume the same number of strange quarks 
and antiquarks in the QGP before hadronization and in the HG after hadronization. This 
number is not expected to be further changed in the HG with local thermal equilibrium. The 
typical thermal hadron momenta are too small for strange hadron production. We expect, 
therefore, the same number of W + Ng in the HG at the chemical freezeout. The total 
entropy is also expected to be conserved approximately during the hadronization of QGP 
and HG expansion. This suggests that the value of Rs at the HG chemical freezeout should 
be close to that in the QGP. The value of Rg for different HG model hts are 0.045 (ID), 
0.037 (VDW), 0.029 (SS) and 0.026 (SS+VDW). It is remarkable that our best HG model 
£t, SS+VDW, where both pion enhancement (r,r = 0, r > 0) and strangeness suppression 
{'js < 1) are included, leads to the value of Rg very close, indeed, to the QGP estimate (P^. 

V. SUMMARY 

In summary, we analyzed new Pb+Pb data at GERN SPS on hadron multiplicities in 
a full dvr-geometry and on antiparticle/particle ratios in the central rapidity region within 
various HG models. The parameters T and ps are found from htting the data and appear 
to be the same, T =165 MeV and /is=235 MeV, within four different versions of the HG 
model. The ideal HG model fails to give a reasonable explanation to the Pb+Pb data sets. 
Meanwhile, the obtained particle number density and total energy density of the ideal HG 
at the chemical freezeout is extremely large, which is inconsistent with the picture of a gas 
of point-like non-interacting hadrons. To improve this situation, the hadron “hard-core” 
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repulsion had been proposed. It removes the artihcially high energy and particle number 
densities in the gas of non-interacting hadrons. The VDW HG model with = 0 and 
Tj = r > 0 leads to the pion enhancement and better £t of hadron multiplicities is obtained. 
Chemical non-equilibrium effects with strangeness suppression were also studied. We have 
found that this modification significantly improves the agreement with data. Furthermore, 
the combined effect of these two mechanisms results in an extremely good agreement with 
the data analyzed. We present the SS-I-VDW HG model predictions for relevant hadron 
multiplicities in central Pb-|-Pb at the SPS. Finally, an interpretation of the obtained results 
in terms of the possible QGP formation at the early stage of A+A collisions at the SPS is 
also discussed. 
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TABLES 

TABLE 1. Particle multiplicities and ratios for central Pb+Pb SPS collisions in different HG 
models. 



ID 

VDW 

SS 

SS+VDW 

Data 

Reference 

T [MeV] 

165 

165 

165 

165 





[MeV] 

235 

235 

235 

235 





V [fm^j 

2490 

6060 

3150 

6640 





r [fm] 

0.00 

0.46 

0.00 

0.45 





7s 

1.00 

1.00 

0.55 

0.62 





Nb 

347 

337 

382 

374 

372+10 


1 


h- 

514 

615 

566 

674 

680+50 

0 

RO 

89.3 

86.4 

60.5 

65.3 

68+10 

0 


11.7 

11.3 

4.46 

5.42 

5.4+0.7 


1 


P - P 

124 

120 

155 

149 

155+20 

0 

K+/K- 

1.78 

1.78 

1.74 

1.75 

1.8+0.1 


E 


p/p 

0.060 

0.060 

0.059 

0.060 

0.07+0.01 

0 


0.242 

0.244 

0.223 

0.227 

0.249+0.019 


E 


Ti/n 

0.496 

0.502 

0.441 

0.452 

0.383+0.081 

1 

21 

1 

xVdof 

108/6 

95.4/5 

12.5/5 

3.58/4 






TABLE II. Particle multiplicities for central Pb+Pb SPS collisions predicted in SS+VDW HG 
model. The values of model parameters are taken from Table I. 


IT 

K+ 

K- 

A 

A 


2 

n 

n 

1850 

83.1 

47.6 

38.9 

4.43 

5.86 

1.33 

0.242 

0.110 
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